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In the limit of infinite spatial dimensions a thermodynamically consis- 
tent theory of the strongly correlated electron systems, which is valid for 
arbitrary value of the Coulombic interaction (U < oo), is built. For the 
Hubbard model the total auxiliary single-site problem exactly splits into 
four subspaces, which describe Fermi and non-Fermi liquid components. 
Such analytical approach allows to construct different thermodynamically 
consistent approximations: alloy-analogy approximation, Hartree-Fock ap- 
proximation and further, that describes the self-consistent renormalization 
of the bosonic excitations (magnons and doublons). 
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In the last decade a lot of the rigorous results in the theory of strongly 
correlated electron systems are connected with the development of the Dy- 
namical Mean-Field Theory (DMFT) proposed by Metzner and Vollhardt 
[1] for the Hubbard model (see also [2] and references therein). It maps 
lattice problem on the effective single impurity Anderson model with the 
generalized partition function 

p = e-/5^orexp l^j^dr f^dr'Y^^r - r')at(r)a.(r')| , (1) 

where C(t{t ~ t') is an auxiliary Kadanoff-Baym field (single-site hopping) 
and for the Hubbard model 

Ho = Un^ni-fi{n'^+ni)-h{n^-ni) = Un^ni-^fi^n^ = ^ ApX^^*, (2) 



* Presented at the XII School of Modern Physics on Phase Transitions and Critical 
Phenomena, L§,dek Zdroj, Poland, June 21-24, 2001 



(3415) 



3416 



A.M. Shvaika 



and there are no restrictions on the U value within this theory. Moreover, 
some classes of binary-alloy-type models [e.g., the Falicov-Kimball model) 
can be studied exactly within DMFT [3]. But in the case of the Hubbard 
model, the treatment of the effective single impurity Anderson model is very 
complicated and mainly computer simulations are used, which calls for the 
development of the analytical approaches [4]. For the partition function (1) 
such approach can be build within a perturbation theory expansion in terms 
of the electron hopping using a diagrammatic technique for Hubbard oper- 
ators [5] , which is based on the corresponding Wick theorem [6] . In present 
article we supplement such strong-coupling approach by the consideration 
of the bosonic fluctuations. 

Consecutive pairing of all off-diagonal Hubbard operators is per- 
formed until we get the product of the diagonal operators only. For the 
single impurity problem (1) and (2) all diagonal operators act at the same 
site, their products can be reduced to the single operators and for the Hub- 
bard model our problem exactly splits into four subspaces with "vacuum 
states" \p) = |0), |2), | |) and | j) and only excitations, as fermionic, as 
bosonic, around these "vacuum states" are allowed. Finally, for the grand 
canonical potential we get [5] 

P p 

where are the "grand canonical potentials" for the subspaces. Now we 
can find the single-electron Green functions by 

Gi^' (r - r') = J^^^^ = E " (4) 

where G<7(p) (''" ~ '^0 axe single-electron Green functions for the subspaces 
characterized by the "statistical weights" 

We can introduce irreducible parts of Green's functions in subspaces S^.^) 

by 

'^;i)M = s;{p)K)-CaK) (6) 

and self-energies in subspaces "^aip) (Dyson equation for the irreducible 
parts) 
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where n^^^^--^ = = for p = 0, cr and 1 for p = 2, o". Here, sclf-encrgy 

^a(j>){^i') depends on the hopping integral Qa'i'^u') only through quantities 

*(t' (p) (<^i/' ) = G^, (p) {iOj,, ) - 'E^i (p) {iOi,, ). (8) 

Now, one can reconstruct expressions for the grand canonical potentials 
r2(p) in subspaces from the known structure of Green functions: 

^(p) = \ - 13 111 [l - Ca('^i/)S(,(p)(w^)] (9) 
P ua 

where ^>(p) is some functional, such that its functional derivative with respect 
to * produces the self-energy: 

l^ j^^t, X =^a(p)K)- (10) 

The second term in (9) corresponds to the sum of the fermionic single loop 

contributions whereas the next ones describe different scattering processes. 
From the grand canonical potential (3) and (9) we get for mean values 

i^a = -^ = ^Wpna(p), (11) 

(0) 1 -r /X ^^(p) 

^^(P) = <(p) + ^a{p){^u) - -g^, 

where in the last term the partial derivative is taken over /x^ not in quantities 
^(8). 

Falicov- Kimball model correspond to the case of Cii'^u) = and, as a 
result, = and S-|-(p)(a;,/) = which immediately gives results of [3] 
(see also [7]). 

Presented above equations allows one to construct different thermody- 
namically consistent approximations. 

The simplest alloy-analogy approximation, which is a zero-order approx- 
imation within the considered approach, is to put $(-p) = 0, T,^^-j(u;i,) = 
and for the Green's function for the single impurity problem one can obtain 
a two-pole expression (Xpg = Xp — Xq) 
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Strong coupling Hartree-Fock approximation takes into account the first 
corrections into the self energy in the form 



^a{p)i<^i^) = U (n^(p) - n^'JJ,)) , (13) 



which gives for the Green's function a four-pole expression: 

(a;,,) = V ^ . (14) 

Now, grand canonical potentials in the subspaces (9) are determined by 
functional 

^(p) = U (n^(p) - n^°^)) (n^(j,) - 4°^^,)) . (15) 

Expression (14) possesses the correct Hartree-Fock limit for the small Co- 

ulombic interaction U and electron or hole concentrations. On the other 
hand, in the same way as an alloy-analogy solution, it describes the metal- 
insulator transition with the change of U. In [5, 8] it was shown that the 
main contributions come from the subspaces p = or 2, that describe the 
Fermi-liquid component, for the low electron (n <C 1) or hole (2 — n, <C 1) 
concentrations, respectively, ("overdoped regime" ofhigh-Tc's) and^)=t,j, 
that describe the non- Fermi- liquid component, close to half filling (n ~ 1) 
("undcrdoped regime"). At low temperatures the Fermi-liquid component 
is in the ferromagnetic (Stoner) state, while the non-Fermi-liquid one is 
antiferromagnetic close to half-filling. 

In order to go beyond the Hartree-Fock approximation one have to con- 
sider, besides the fermionic excitations, also the bosonic ones [5] which cor- 
respond to the creation and annihilation of the doublons (pairs of elec- 
trons), described by the X^'^ and X^^ operators, for subspaces p = 0,2 and 
magnons, described by the X^^ and X^'^ operators, for p =1, j. The single 
loop contributions of such bosonic excitations can be summed up and one 
can obtain 

^ip) = ^^^^[^ ~ ^ ^ ^ ^<Ta(p) ('^m)) Xaaip) ('^m) , (16) 



where 



Daa(p){i^m) = f ; ^ "T, (17) 



Xaaip){^m) = --^ ^ ^a(p) ('^m-i/) (18) 
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for subspaces p = 0,2 and 

Daa{p){^m) = : p ^ ^T, (19) 

Xaa{j,){^m) = ~Y^-^ a{p){'^v-m) (20) 



for subspaces p = a,a, which gives for mean values the following expression 
= + (21) 



P m 1-U(l ±UD^^^p){uJm)) Xaa{p){^m) 

Now self-energy contains the frequency dependent part 

^a(j>) = U (n^(p) - nfl^) + S,(p) (a;,), (22) 

ol ^ (1 ± t^-Do-CT(p)(<^ 

= ^ 75Z^^ 7 ^ X 

P m l-U \l± UD^^f^p) (uJrn) ) Xaa{p){'^m) 

X / *<^(p)(^m-i.), for P = 0, 2 

I ^a{p)ii^u-m), for p = cr, ' ^ 

that describes the contributions from the doublons for the Fermi liquid com- 
ponent {p = 0, 2) and magnons for the non-Fermi liquid one {p =t, J,) with 
the renormalized spectrum determined by the zeros of denominator in (23). 

Expression (16) for functional has the same form as the correction 
to free energy in the theory of the Self- Consistent Renormalization (SCR) 
of spin fluctuations by Moriya [9]. But in our case it describes contributions 
from the single-site bosonic (magnon or doublon) fluctuations with specific 
renormalization functions different for different subspaces. Spin fluctuations 
give the main contribution close to half filling in the non- Fermi liquid regime 
but for small electron (n <C 1) or hole (2 — n <C 1) concentrations the 
contributions from the doublon (charge) fluctuations must be taken into 
account. 
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